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Abstract. The energy spectrum of the electromagnetic normal modes plays 
a central role in the theory of the van der Waals and Casimir interaction. In 
dissipative media it appears as if there are distinct normal modes with complex 
valued energies. The summation of the zero-point energies of these modes render 
a complex valued result. Using the contour integration, resulting from the use of 
the generalized argument principle, gives a real valued and different result. We 
resolve this contradiction and show that the spectrum of true normal modes form 
a continuum with real frequencies. We illustrate the problem in connection with 
the van der Waals interaction between two metal half spaces. 
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1. Introduction 

The interaction energy in a system may be expressed in terms of the energy shifts of 
the electromagnetic normal modes of the system [I] ; the longitudinal bulk modes can 
be used to find, e.g., the polaron energy in a polar semiconductor or the exchange- 
correlation energy in a metal; the transverse modes to find the Lamb shift [2]; the 
surface modes to find the van der Waals interaction between objects; the vacuum 
modes to find the Casimir [3] interaction. In all these cases the energy, or frequency, 
of a normal mode is found as a solution to an equation of the type 



the condition for having a mode. The function / (u), which involves the dielectric 
properties of the system, is often obtained as a determinant of a matrix. The solution 
to equation (fT]) approaches a pole of the function / (ui) when the interaction in the 
system is gradually turned off. When the interaction is turned on the zeros move 
away from the poles and fi/2 times this shift in frequency is the contribution to the 
interaction energy of this particular mode; the interaction energy is the change in total 
zero-point energy when the interaction in the system is turned on. 

It is often straight forward to use this sum-over-modes approach, or mode- 
summation method, to obtain the interaction energy as the following sum over zeros 
and poles of the function / (w): 



/ M = o, 



(i) 




(2) 
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However, sometimes it is not possible to find an explicit solution to equation ([T]) and 
sometimes the poles and zeros form a continuum. In those cases the result may be 
obtained with the so called generalized argument principle pQ, 



where the integration is performed along a contour including the poles and zeros in 
the right half of the complex frequency plane. The integration should be performed 
in the positive sense, i.e. in the counter-clockwise direction. 

The specific problem of the force between two semi-infinite dielectric slabs 
separated by a dielectric layer has been solved with different approaches; in Lifshitz' 
[3] very complicated theory the dielectrics are characterized by randomly fluctuating 
sources as demanded by the fluctuation-dissipation theorem; Schwinger et al. [5] 
derived the force using Schwinger 's source theory, where "the vacuum is regarded as 
truly a state with all physical properties equal to zero" ; van Kampcn et al. [6] applied 
equations ((2|) and ((3]) to the surface modes of the geometry to find the interaction 
energy and force. 

All these approaches, that appear to be quite different, give one and the same 
result. Milonni and Shih [7j made a study, on how these theories are related, based 
on conventional QED. The result in equation ([3]) is consistent with what one arrives 
at from many-particle theory where the focus is put on the interacting particles in 
the system, not the electromagnetic normal modes; there is no explicit reference to 
zero-point energies. This is demonstrated in detail in pQ in the case of the exchange- 
correlation energy of a metal. In [8] the van der Waals and Casimir forces between 
two quantum wells were derived both in terms of the zero-point energy of the normal 
modes and as the result from correlation energy; both approaches produced the same 
result. Thus, there are different complementary approaches to the interactions in a 
system. 

In the present work we are concerned with the sum-over-modes approach in 
presence of dissipation in the system. Van Kampen et al. [6] considered non- 
dissipative materials, only, with real-valued dielectric functions. In the case of 
dissipative materials the sum-over-modes approach runs into yet another problem. 
The straight forward solution of equation (TT]) produces complex- valued frequencies or 
energies, signalling that the normal modes are no longer stable - the modes decay. If 
this were the case they would no longer be true normal modes. The interaction energy 
in equation ((2]) becomes complex valued. Using instead equation (J3j) produces a real- 
valued interaction energy. This situation has made many researchers confused and 
led to believe that there is something fundamentally wrong with the mode-summation 
method. With the present work we will try to resolve this confusion and show that the 
result of equation ([3]) is the correct result in the case of dissipation; the complex-valued 
result of equation (|2|) is not correct. In an earlier brief report [9] we have illustrated 
the problem by considering the longitudinal bulk modes in a metal. In the present 
work we study the modes associated with the van der Waals interaction between two 
metal half spaces. 

In section[2]we discuss the analytical properties of the dielectric function and show 
that its zeros and poles are all located at the real frequency axis. In section[3]we apply 
the different formalisms to the specific example viz. the van der Waals interaction 
between two metal half spaces. We introduce a modified mode-summation method 
that works reasonably well in section 0] and finish with summary and conclusions in 
section [5j 




(3) 
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2. Analytical properties of a dielectric function 

In our illustrating example we will use metals and represent the dielectric function 
with one of the Drude type, 

e(w) = l-w2,/w(w + t»7). (4) 

We let the parameter rj be a positive real-valued constant; n, which is the result 
of electron scattering against impurities or other defects, is really w-dependent and 
complex valued; the approximation we use is good for a metal in the low momentum 
limit, for frequencies below the plasma frequency, ui p i. 

Before we continue let us discuss the general analytical properties of a dielectric 
function. The physical dielectric function, the one that can be measured in 
experiments, exists on the real frequency axis, only. It is retarded, which means that 
it obeys causality. In theoretical treatments one obtains a function that is analytical 
in the whole complex frequency plane except on the real axis, where all the poles are 
situated. To obtain the retarded version one either shifts all the poles downwards to 
an infinitesimal distance below the real axis and perform the calculation on the real 
axis. Alternatively one lets the poles stay put and perform the calculation just above 
the axis. For the discussion in this work it is better to use this last method. There 
are other versions of the function, advanced, timeordered and antitimeordered. The 
timeordered is often used in many-body calculations since it allows some very useful 
theorems to be used. With this version one calculates the function just above the 
positive real axis and below the negative real axis. All different versions are identical 
everywhere except at the real frequency axis. From now on, if not stated otherwise, 
when we discuss the dielectric function we mean the function with its poles on the 
real axis. The function has the properties s (— w) = e (lu) and e {to*) = e {to) ■ From 
these equations follows that the relation between the different forms of the function 
in the lower and upper half plane is ei (lu) — e u (u)*) . For the full dielectric function 
with frequency dependent rj these two analytical expressions are the same but not 
when -q is treated as a constant. Then we have e u (lu) = 1 — to^j lu (lu + irj) and 

ei (lu) = 1 — Lu 2 pl I lu (lu — irf), respectively. In the full treatment we have 

irj (lu*) — [irj (lu)]* = —i [rj (lu)}* = — iRc [rj (lu)] — Im [rj (lu)] , (5) 

and we may identify 

Re [rj (lu*)] = -Re [77 (lu)] ; Im [7/ (lu*)] = Im [7? (lu)] . (6) 

We see that the dominating real part changes sign when we cross the real axis while 
the imaginary part does not. 

As an illustration, let us look at the expression for the dielectric function of an 
impure metal in the so called generalized Drude approach [10) . Let nj, n, S(q), 
and LUo(q) be the density of impurities, electron density structure factor for the 
impurities and impurity potential, respectively Then the dielectric function in the 
small momentum limit is 

e(iu) = l-— ^ p (7) 

" I " + afei E S (q) <Z 4 N (q) I 2 [e- 1 (q, u>) - ^ (q, 0)] J 

where f2 is volume of the system. The second term within the parentheses in the 
denominator is irj. For a large system the summation over momentum is usually 
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replaced by an integral over a continuous momentum variable. Let us now keep the 
discrete summation. The derivation of the dielectric function of the pure metal at finite 
momentum also contains a discrete summation, now over the electron momentum. The 
function in RPA (Random Phase Approximation) is 

e(q, w) = 1 + («,/n) E n(k) [1 - n (k + q)] 

x + (e k+q - eiOP 1 - [nu - (ek+q - £k)] _1 } , 

where n(k) is the Fermi-Dirac occupation number. We see that, if we make the 
calculation just off the real frequency axis, the imaginary part consists of a sum of 8- 
functions, infinitesimally spaced when the volume of the system goes to infinity. These 
form the single particle continuum in the cjq-plane. The real part passes through zero 
between each neighbouring pairs of 5-functions. When the volume goes to infinity 
one can replace the summation by an integral. The imaginary part then turns into a 
smooth continuous function and the real part does no longer pass through zero inside 
the continuum. When one wants to find the zeros and poles of the function e(w), in 
equation ([7]), one should keep the discrete summations everywhere. Then one realizes 
that this function, also, has its poles and zeros on the real axis and that they are in 
the form of a continuum. 

If we now instead let all summations turn into integrals the parameter 77 is a 
smooth complex valued function of frequency. However its real part dominates and is 
almost constant for frequencies below the plasma frequency. From the expression in 
equation ([7]) we can easily verify that the real part of rj changes sign when one crosses 
the real axis. For higher frequencies the contribution to the interaction energy quickly 
drops off so equation ((4]) is good enough for our purpose here. 



3. An illustrating example 

As an illustration we have chosen the van der Waals interaction between two metal 
half spaces separated by the distance d. The modes in this system are characterized 
by the 2D (two-dimensional) wave vector k. There are several modes for each wave 
vector. To find the modes of wave vector k we let f{oS) in equation (HJ be [T] 

/ k M = [e M + l] 2 - e- 2W [e [u>) - l] 2 = {[e (u>) + 1] - e~ kd [e ( w ) - 1]} 

x { [e (a,) + 1] + e- kd [e ( W ) - 1] } = ft (u) ft (co) . (9) 
Equation ([T]) then gives us four zeros, The two in the right half plane are 

( u pl J[l + coth (kd/2)}- 1 - ( V /2cu pl ) 2 - t V /2; ft (w) = 

^k=< \ (10) 

[ w p/v /[l + tanh(fcd/2)] 1 -(ri/2u> pl y -irj/2; ft(u) = 

Thus the zeros are below the real frequency axis. The problem is that the expression 
in equation ((4} for the dielectric function is only valid above the real axis. Below the 
real axis 77 has the opposite sign. The zero is very illusive. If one approaches the zero 
from the upper half plane it makes a jump to the upper half plane when one crosses 
the real axis. This clearly shows that it is not safe to use equation @, directly. If 
plotted the real part of the mode energies from equation (TTUj) form two branches of 
modes, one acoustical (lower) and one optical (upper). The lower (upper) branch 
comes from putting ft = (ft = 0). A first naive guess is that a good approximation 
of the interaction energy would be to just sum the real part of the zero-point energies 
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for the modes. This might seem to be a good and simple short cut to an approximate 
result. However, as we will see later, it turns out not to be such a good idea. 

Let us now instead use equation to find the energy contribution from each 2D 
wave vector. The total interaction energy per unit area is 



E 



ft ^ 

k 



E k 



-E- 

ft 4^ 2ni 



dz—z — In fk (z), 
2 dz J V ; ' 



(11) 



where ft is the total area of one half space. In an actual calculation we let this area 
go to infinity and the summation over k turns into an integration. We first choose our 
contour to encircle the positive real axis. The contour then consists of two parts; an 
integration from zero to plus infinity performed just below the real axis; an integration 
from plus infinity to zero performed just above the axis. We may first perform an 
integration by parts in both contributions and end up with 



1 Tb f 

Then the integration below and above the axis are combined into one integral 



1 h f du fk(coy 



E u = — - / —In 



i2j2TT / k (w) 



2 

E 

1=1 



OO 

\ J ^tan-'ilmlfU^/ReifU")]}, 



(12) 



(13) 



where the function tan^ 1 is taken from the branch where < tan -1 < ir. 

Let us now instead deform the integration contour into a semicircle in the right 
half plane with the center of the circle at the origin and the straight part parallel with 
and just to the right of the imaginary frequency axis, 



E k 



- / — 21 

2ir 



= - / — 21 

2ir 



(14) 



The integration along the curved part of the contour vanishes when the radius goes 
to infinity. The four in the denominator takes care of this. The zeros and poles 
of the function /k (w) in equation §§§ do not change from the division by a factor 
of four. Thus we end up with two integrals for the same thing; one along the real 
frequency axis; one along the imaginary axis. Now, we may simplify the integrands by 
letting x = (ri/2)/ujpi and expressing the frequency in terms of the plasma frequency. 
Equations p^|) and (fLlj) then reduce into 



E v 



huo p i/2 



OO 


tan 1 














2x (1 



uj {2lu 2 + 8x 2 - 1 



+ tan 



2x1 



to (2uj 2 + 8x 2 - 1 - e~ kd ) 



dtoF (uj), 



and 

gk 

hujpi/2 



duj nl \2 + 1/uj(uj + 2x)V 

— 2ln 

2tt 



e- 2kd [1/lu (lo + 2x)Y 



dioG (uj) 



(15) 



(16) 



respectively. These two integrands F (uj) and G(u>) are shown in figure [T] for a wave 
vector with kd = 0.5, for the two parameter choices x = 0.1 and 0.01. 
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Figure 2. The contribution to the interaction energy from one wave vector with 
kd = 0.5. The solid curve is the exact result as calculated from either equation 
(1 1 5 I t or II16I I . The dashed curve is from neglecting the imaginary parts of the poles 
and zeros. The circles are the result from using equation I118H and the parameters 
given in the text. 



In figure we present the results from the different approaches. The solid curve 
is the exact result from equation (fT5|) or (| 16|) and the dashed curve shows the short 
cut from just dropping the imaginary parts of the mode energies. The circles are the 
result from yet another approximation that we will discuss in section 2) The fact that 
equations (fT5)l and (|16[) produce the same result proves that the poles and zeros are 
all on the real axis. 
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4. A mode summation method that works 



In the preceding section we found that the naive use of the mode-summation method 
failed. The question is whether the mode-summation method can be used at all in the 
case of dissipative materials. To investigate this we make a Lehman representation of 
the dielectric function, 



e(w) = 1 + 



dbJ 



— 2Ime ( uj 



2tt 



= 1 - 



dio 



2tt 



{uj/ujpi) - (a/) (u/) +4x 2 



(17) 



The integral over frequency can be viewed as a limit of discrete frequency summations 
where the step size goes towards zero. In each summation the dielectric function has 
its poles and zeros on the axis. These points come closer and closer when we take 
the limit. The integrations over frequency and momentum should always be viewed 
as limits of discrete summations since the system is always finite in size and possible 
energy and momentum transfers are discrete. 

Now if we approximate the integral of equation (|17[) with a discrete summation, 



e(w)wl 



8x LJ„ 



2tt i n 



£ (^M>0 2 -(^) 2 (^) 2 + 4x 5 



i=l 



(18) 



where u) i — (i — 1/2) uo max j i max , the dielectric function has a finite number of poles 

and zeros. Using this function in the expressions for the two factors (uj) and (u>) 
of equation © means that also these factors each have the same number of poles and 
zeros. Using equation |(5J) where the sum runs over these poles and zeros gives us an 
approximate result for E^, the contribution to the interaction energy from mode k. 
The result becomes asymptotically the exact result when we let w max and i max go 
towards infinity. In figure [2] the first (last) two circles were calculated including 50 
equidistant poles in the region below 2w p i (3w p i), i.e. u) max = 2(3) and i max = 50(50). 
Thus we have to find the frequency of 50 poles and 50 zeros. This is feasible and the 
result is in much better agreement with the exact result than what one obtains when 
just neglecting the imaginary part of the zeros and poles. 

In this treatment all poles and zeros are real-valued. We now have one acoustical 
and one optical branch in between each pair of neighbouring poles. When the 
interaction is turned off the zeros end up at the position of the lowest of the two 
poles. When the interaction is turned on they move up a distance and it is these 
shifts in frequency or energy that are their contributions to the interaction energy. In 
figure [3] we show as open circles these shifts for the acoustical (the curve with left- 
most maximum) and optical branches. The combined shifts are shown as filled circles. 
We clearly see the resemblance of the filled circles to the integrand of equation JTJ 
represented by the solid curve with circles in figure [T] 



5. Summary and conclusions 

In summary, we have demonstrated that the energy of the electromagnetic normal 
modes in dissipative media are real valued and lead to real- valued interaction energies 
and forces. The modes appear to be distinct with complex valued energies, but they 
are not. The modes form a continuum on the real frequency axis and the direct use 
of the mode summation method is no longer feasible. The functions appearing in 
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Figure 3. The energy shifts of the zeros away from the poles when the interaction 
is turned on. The open circles are the shifts of the modes corresponding to one 
of the two branch types. The filled circles is the combined effect. Sec the text for 
more details 



the condition for modes contain integrals over momentum or frequencies. Since the 
system is finite in size these integrals should be considered as discrete summations; the 
possible momentum and energy transfers in a finite system is discrete. In doing so all 
zeros and poles end up on the real axis. When one takes the limit when the volume goes 
to infinity the zeros and poles form a continuum of points on the real axis. We have 
demonstrated two different contour integrations that can be used to find the result; 
the one we favour most ends up with an integration along the imaginary frequency 
axis; the other involves an integration along the real frequency axis. We have further 
demonstrated that one may use an approximation in which this continuum of modes 
is replaced by a finite number of distinct modes and obtain reasonable good results 
using the mode-summation method. This gives a great improvement compared to the 
result from just neglecting the imaginary parts of the complex valued energies of the 
apparent modes. 
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